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Abstract 

As helium-4 is cooled below 2.17 K it undergoes a phase transition to a fundamentally new state 
of matter known as a superfluid which supports flow without viscosity. This type of dissipationless 
transport can be observed by forcing helium to travel through a narrow constriction that the normal 
liquid could not penetrate. Recent experiments have highlighted the feasibility of fabricating 
smooth pores with nanometer radii, that approach the truly one dimensional limit where it is 
believed that a system of bosons (like helium-4) may have startlingly different behavior than in three 
dimensions. The one dimensional system is predicted to have a linear hydrodynamic description 
known as Luttinger liquid theory, where no type of long range order can be sustained. In the 
limit where the pore radius is small, Luttinger liquid theory would predict that helium inside the 
channel behaves as a sort of quasi-supersolid with all correlations decaying as power-law functions 
of distance at zero temperature. We have performed large scale quantum Monte Carlo simulations 
of helium-4 inside nanopores of varying radii at low temperature with realistic helium-helium and 
helium-pore interactions. The results indicate that helium inside the nanopore forms concentric 
cylindrical shells surrounding a core that can be fully described via Luttinger liquid theory and 
provides insights into the exciting possibility of the experimental detection of this intriguing low 
dimensional state of matter. 



I. INTRODUCTION 



As the spatial dimension of an interacting many-body system is reduced, both thermal 
and quantum fluctuations are enhanced, providing a fascinating arena for the study of 
complex phenomena. It is well known that in strictly one dimension, there is no broken 
continuous symmetry but instead the persistence of only quasi-long range order characterized 
by power-law decay of correlation functions. The simplest example of a model displaying 
these features is the one dimensional (Id) non-relativistic Galilean invariant Bose gas with 
delta function interactions. The ground state of this model has been known for almost fifty 
years and can be solved exactly using the Bethe AnsatzR An alternative approach, based 
on the concept of a low dimensional harmonic fluicP^ was first understood to be universal 
by Haldane^ with the techniques formalized in these seminal works now generally known 
as bosonization and Luttinger liquid (LL) theory. The combination of these two approaches 
indicate that there is no well defined quantum phase transition at zero temperature as a 
function of the strength of the interactions in Id. As interactions are increased, the system 
crosses over from a regime dominated by phase fluctuations, to one with tendencies towards 
density wave order. Likewise, there is no phase transition as a function of temperature. 

Historically, the experimental study of low dimensional physics has been confined to 
fermionic systems such as insulating spin chains and quantum wires^ and a bosonic real- 
ization has been lacking. The identification of accessible bosonic systems is particularly 
relevant in light of the universal predictions of LL theory. Recently, weakly interacting 
ultra-cold bosonic gases consisting of laser trapped Alkali atoms confined to cigar shaped 
quasi-one-dimensional condensates^ or rings^ have been studied but it is difficult to obtain 
results at higher densities where the constituent atoms are more strongly interacting. The 
various theoretical approaches and considerations required to study such systems have been 
covered in detail in a recent review article^. 

A considerably different approach involves the preparation of a quantum fluid of interact- 
ing bosons, such as neutral 4 He in a physically confining geometry at high density. In fact, 
the ability of helium-4 cooled below T\ = 2.17 K (the bulk superfluid transition tempera- 
ture in 4 He) to flow through a narrow constriction (superleak) is one of the original defining 
characteristics of the superfluid state of matter. Original experiments in this vein achieved 
physical confinement of helium through the quasi-lD cavernous networks of porous glasses 
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such as Vycor^l and more recently using folded sheets of mesoporous materials^. Advances 
in nanofabrication techniques have culminated in a more systematic approach that has been 
exploited by Savard et a/ J 14 * 15 * in studying 4 He inside nanopores of variable radii by sculpting 
a pore through a SisN4 membrane using a transmission electron beam. Although these ex- 
periments have thus far focused on flow properties of helium in the gas phase in nanopores^ 
and the superfluid phase for wide nanoholes^ they provide a tantalizing road map for the 
experimental detection of a bosonic Luttinger liquid. It is thus natural to ponder a simpler 
equilibrium system, that can be studied numerically as a function of pore radius and tem- 
perature below T\. We expect that when the pore radius is sufficiently small, the length 
sufficiently long and the temperature low enough with respect to T\ the 4 He filled nanopore 
should begin to exhibit strictly Id behavior where LL theory can be used to characterize 
the nature of the quasi-long-range order. 

The remainder of this paper is concerned with the precise definition of sufficiently in 
the previous sentence and is organized as follows. We first identify the particulars of a 
theoretical model of confined helium-4 and introduce the stochastically exact numerical 
method employed to study its behavior. The numerical results are used to determine a phase 
diagram for helium inside the pore with special attention payed to the types of structures that 
are allowed by the Galilean invariant cylindrically symmetric confining potential. Density- 
density correlations inside the inner region of the pore are then analyzed in terms of Luttinger 
liquid theory which is introduced in an economical fashion with references provided to more 
complete and elaborate treatments elsewhere. The results highlight the applicability of the 
linear hydrodynamics of LL theory in describing helium inside nanopores and a discussion 
of the consequences and limitations of this finding is presented. 

II. PATH INTEGRAL SIMULATIONS OF CONFINED HELIUM-4 

The finite extent and translational invariance of the van der Walls interaction between 
4 He atoms combined with their fundamental indistinguishably conspire to make the numer- 
ical study of quantum fluids a challenging task. Unbiased stochastically exact simulations 
of the quantum bosonic many-body problem in the continuum at low temperature have 
only recently become more feasible through the introduction of the continuous space Worm 
Algorithm (WAP 1121 . This method builds upon the conventional Path Integral Monte Carlo 
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approach of Ceperley^ that exploits the Feynman path integral interpretation of quantum 
mechanics to perform Metropolis sampling of the worldlines of particles in d + 1 dimensions. 
The WA extends the configuration space to include open worldlines known as worms that 
are not periodic in imaginary time and directly contribute to the single particle Matsubara 
Green function allowing for simulations to be performed in the grand canonical ensemble. 
The simulation cell is kept in thermal equilibrium with a surrounding bath with which it 
can exchange particles and the chemical potential \i is an input parameter providing access 
to important physical observables such as the compressibility. In addition, the WA employs 
an efficient method via the swap operator^ to directly sample the quantum statistics of 
identical bosons through only local exchanges of worldlines. Such exchanges are the crucial 
step en route to the superfluid state of matter, with the superfluid density being measured 
through a topological winding number; the standard technique in continuum simulations^. 

We wish to perform WA simulations on a realistic physical model of the experimental 
single nanopore geometry discussed in the introduction with details in References HU and 
T51 and thus must consider interactions between helium atoms, as well as the interactions 
between helium and the atoms composing the material surrounding the nanopore. Our 
starting point is the general microscopic many-body Hamiltonian 

- r i i 



i<3 

where are the spatial positions of the N helium atoms and we have set h — k# — 1. The 
natural units of length in the nanopore systems are angstroms and in these reduced units, 
the mass of the helium atoms is m = 0.0826 A K~ x . The interaction energy between two 
helium atoms v, is given by the 1979 Aziz potential 

v(r) = tu (—] , (2) 

<\ 4— (Ce.Cs C 10 \ Je-W- 1 )' ; x<D 
u{x) = Ae -[— + — + — )< (3) 
\x e x 8 x 10 J | 1 ; x>D 

with a set of parameters inferred from both first principles calculations and experimental 
measurements to be r min = 2.9673 A, A = 5.449 x 10 5 , e/k B = 10.8 K, a = 13.353, 
D = 1.2413, C 6 = 1.3732, C 8 = 0.42538 and C 10 = 0.1781. This potential, shown as an 
inset in Figure [TJ contains both hardcore repulsion for distances less than r^ c ~ 2.64 A, a 




FIG. 1. The interaction potential V between helium and the walls of the nanopore constructed of 
amorphous S13N4 for various radii ranging from R = 2.9 A (left) to R = 12.0 A (right) calculated 
via Eq. Q where r is the radial distance of a helium atom from the axis of the cylinder. The inset 
shows the Aziz interaction potential v, described in Eq. ^ for two helium atoms separated by a 
distance r. 

minimum at r min and weak attraction at large separations. The external potential, V can 
be computed by modeling the nanopore as long cylindrical cavity carved inside a continuous 
medium. The potential energy of interaction of a single helium atom with the atoms of 
the medium can be obtained by integrating a Lennard- Jones pair potential as described in 
Reference [21] yielding 

i „. TTE(7 3 n 17 /cr\ 9 rr\ /cr\ 3 /r\ 



(4) 



for a perfect cylindrical cavity of radius R where n is the number density and a and e are 
the Lennard- Jones parameters of the surrounding medium. The functional coefficients are 
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given by 

ug(x) = 1 — [(1091 + 11516a; 2 + 16434a; 4 + 4052a; 6 + 35x 8 )E(x 2 ) 

-8(1 - a; 2 )(l + 7a; 2 )(97 + 134a; 2 + 25a; 4 )ir(a; 2 )] (5) 
u 3 (x) = _ 2 [(7 + x 2 )E(x 2 ) - 4(1 - x 2 )K{x 2 )] (6) 

where K(x) and E(x) are the complete elliptical integrals of the first and second kind re- 
spectively and r is the radial distance from the axis of the pore. Using the Lorentz-Berthelot 
mixing rules for the microscopic Lennard- Jones parameters of amorphous S13NJ22H2I we se ^ 
e = 10.22 K and a = 2.628 A with n = 0.078 A -3 . A plot of the pore potential for different 
radii is shown in Figure [I] The use of parameters for different media such as silica glass^ 
will alter the overall energy scale and may change the excluded volume experienced by atoms 
inside the pore but will not qualitatively affect the findings of this study. 

III. HELIUM-4 INSIDE NANOPORES 

Using this model for the interaction and confinement potentials, we have computed the 
low temperature thermodynamics of fluid 4 He inside nanopores with lengths between L = 50 
and 200 A and temperatures ranging from T = 0.5 — 2.0 K using the WA. We imagine the 
nanopore to be immersed in an essentially infinite bath of helium maintained at saturated 
vapor pressure (SVP) which sets the chemical potential in our grand canonical simulations 
to be fi = —7.2 K. A wide range of cylindrical pores with radii R between 2.9 and 12.0 A 
have been considered where we assume periodic boundary conditions along the axis of the 
cylinder. 

We restrict the temperature in our simulations to be less than T\ ~ 2.17 K and wish to 
confirm that we are studying a low energy quantum regime. This is most straightforwardly 
determined by measuring the average kinetic energy per particle in our quantum Monte 
Carlo (QMC) simulations as shown in Figure [2] and comparing with the temperature. We 
observe a general trend of decreasing kinetic energy, from more than 30 K per particle for 
R = 2.9 A to 20 K per particle for R = 12 A. The non-monotonic behavior observed between 
R = 2.9 A and R = 4.0 A can be attributed to the complete transverse confinement of helium 
in the narrowest pore. The minimum value of the kinetic energy observed is an order of 
magnitude larger than all temperatures considered and we thus conclude our simulations are 
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FIG. 2. The average kinetic energy per particle vs. pore radius measured via a thermodynamic 
estimator at low temperature (T = 0.5 K) for a pore of length L = 100 A held at fi = —7.2 K. 

dominated by quantum effects. Unless otherwise specified, we fix the length of the nanopore 
to be L = 100 A allowing three dimensional number densities to be easily converted into 
particle numbers (usually near 1000 helium atoms). The errorbars in Figure [2] are the result 
of a bootstrapping analysis of over 10 6 individual measurements of the kinetic energy using 
a thermodynamic estimator^. We have fixed the number of time slices (the discretization of 
the imaginary time interval corresponding to the inverse temperature) to be 125 per degree 
kelvin. Possible Trotter error introduced by using a finite imaginary time step At has been 
analyzed and we find that it is well described by a term proportional to (At) 4 consistent 
with the use of a fourth order path integral factorization of the density matrix^. 

Assured we are in a low energy regime over the entire range of temperatures considered 
in this study, we now switch our attention to the structures formed by helium atoms inside 
the pore. There exists a large literature on the types of phases that are formed when 
confining a quantum fluid of helium in a cylindrical cavity and a complex phase diagram has 
been predicted for carbon nanotubes^ 3 ^ anc [ both smoottp3H22] and porous nanopore J 36 * 37 l 
Advances owing to the development of the WA have allowed this and a previous worfP^l to 
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study considerably larger radius pores at high density and finite temperature in the grand 
canonical ensemble that approach those studied in recent experiment J 14 * 15 ! 

For fixed chemical potential p = —7.2 K the average volume density py = (N)/itR 2 L 
inside the nanopore is shown in Figure pi where (• • • ) indicates a thermodynamic average 
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FIG. 3. The average volume density as a function of nanopore radius for helium at saturated vapor 
pressure and T = 0.5 K. 

performed in the quantum Monte Carlo. The observed general trend is an increase in the 
volume density as a function of radius continuing to a limiting value for large R that is 
approaching the expected bulk density of 0.02198A The resulting phases of helium 
inside the nanopore can be elucidated by measuring the radial density of particles, Pr{t) 
defined to be the number of particles found a distance r from the axis of a pore of radius R 
normalized such that the linear density pi is 



PL 



(N) 



2it J drrp R {r) 



(7) 



where pn(r) includes an implicit average over the axial and angular coordinates. The results, 
shown in Figure [4] display a progression of structures including a nearly one dimensional 
chain of helium atoms for R = 2.9 A to a series of three concentric shells surrounding a 
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FIG. 4. The radial number density pn(r) of helium inside nanopores with radii R = 
2.9,4.0,6.0,8.0,10.0,12.0 A with an inset showing a full instantaneous worldline configuration 
inside the pore measured in the quantum Monte Carlo simulations projected on z = 0. All pores 
of L = 100 A and are held at saturated vapor pressure corresponding to fi = —7.2 K. 



chain for R = 12.0 A. The insets display a snapshot of the discretized helium worldline 
configuration (space and imaginary time coordinates) from the QMC simulations projected 
into the xy-plane at z = 0. 

The observed oscillations in the radial density can be easily understood by contemplating 
the filling of an empty nanopore through a classical analysis of the potential interactions 



only (Figure [I]). Beginning with the widest pores, we observe an excluded volume effect due 
to the hard Si3N 4 wall, as well as a deep minimum of the potential near the surface. This 
will naturally lead to an adsorption effect or wetting of the pore surface with helium atoms 
forming a shell due to the cylindrical symmetry. As helium atoms continue to enter the pore, 
the adsorption will cause the area density of the shell to increase until the average separation 
between atoms inside the shell begins to approach the hard core radius r^ c of the interaction 
potential. At this point, it will become energetically favorable to form a new layer inside 
the one adsorbed on the surface. As this process continues, a series of concentric cylindrical 
shells may be formed inside the pore analogous to the layering observed in quantum films 
of bosons^. 

The results of the radial density for pores of different radii can be separated into two 
natural groups based on the presence of helium at high linear density near the center of the 
pore. The existence of such an inner cylinder (IC) which can be thought of as a quasi-lD 
chain of helium atoms for a given radius depends on the details of both the helium-pore 
and helium-helium interaction potentials. The pore wall (and thus the details of the wall 
material, SisN 4 here) produces a region of excluded volume, forcing the helium atoms to 
remain a distance greater than r e ~ 1.32 A from the walls of the pore. This sets the location 
of the outermost shell. The separation between shells is then restricted to be near the 
minimum of the Aziz potential r m ~ 3 A with some weak dependence on the radius of the 
pore due to screening of the inner shells. We find that an IC exists whenever the pore radius 
is approximately divisible by three: R = 2.9, 6.0 and 12.0 A. 

The presence of an IC of helium is intriguing, as we expect that such a quasi-ld bosonic 
system should lack any long range order down to zero temperature and be described by 
Luttinger liquid theory. In the next section we test this prediction by introducing the 
universal theory for an effective harmonic fluid and use it to exhaustively study nanopores 
with radii R = 2.9 A and R = 12.0 A. 

IV. A LUTTINGER LIQUID CORE 

In an attempt to understand the relevant low energy degrees of freedom for quasi- ID 
helium atoms inside the pore, we begin by studying a microscopic Hamiltonian like the one 
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in Eq. ([!]) in second quantized form for a strictly one dimensional system of bosons 

H = J dz ^n & ^^ d ^^ + \f dz'p{z)v lA {z-z')p{z) (8) 



where ^(z) = \J p(z)e is a bosonic creation operator with ^^(^r 7 )] = 5(z — z') 

and the Id density and phase operators satisfy 



[p(z),e^ z,) } = 5(z- z')e l ^ z) . 



(9) 



It is the ability to simulate such microscopic Hamiltonians over a range of energy scales that 
makes the WA so attractive. The stochastically exact equilibrium properties of a system of 
strongly interacting bosons at finite temperature can be determined using only fundamental 
input parameters such as the particle mass m and the details of the interaction potential. 



A. Luttinger Liquid Theory 

The manipulations (generally referred to as bosonization) that transform a microscopic 
one dimensional interacting Hamiltonian like Eq. ^ into a universal description of the linear 
hydrodynamics of a quantum fluid (either bosonic or fermionic) are by now standard 6 1 8 1 11 1 39 1 . 
In particular, for bosonic systems, Reference 1391 provides a detailed and pedagogical deriva- 
tion of the Luttinger Liquid Hamiltonian defined by 

HiXj — pN = J dz [ Vj {d z( j>f + v N {d z 9f] , (10) 

where we have included only those terms which are formally relevant in the renormalization 
group sense. The velocities vj and v jv are fixed by the microscopic details of the underlying 
high energy model and the angular field 6(z) appears through the redefinition of the density 
operator 



Po + -d z 9(z) 

71 



J2 e ^me(z) (n) 



where po is the number density at T = in the thermodynamic limit. The underlying 
bosonic symmetry requires that Eq. (11) in conjunction with Eq. (|9| produces the following 



commutation relation 

ld z B{z),cj>{z')]=i>!T5{z-z'). (12) 

11 



If the system described by Eq. (10) exhibits Galilean invariance we can identify 



V N 



m ' 
1 



(13) 

(14) 



where k is the adiabatic compressibility in the limit L — > oo, T — > C 6 . The connection 
of these effective velocities to the underlying microscopic details of the original interact- 
ing Hamiltonian is now clear: a highly incompressible state with vn ^> 1 should have a 
nearly constant density and thus d z 6(z) ~ whereas a state displaying phase coherence has 
d z (j){z) ~ and thus at T = there should be a finite superfluid fraction with vj 3> 1. 

For a system with periodic boundary conditions in the ^-direction, our original boson 
field must satisfy W(z + L) = W(z) leading to a mode expansion of 9(z) and <j)(z) indexed 
by wavevector q = 2~Kn/L where n is an integer® 
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(--(6j + 6_,)sgn(g), (15) 
e^(6t-6_ ? ). (16) 



Here, 6j (6 g ) is a bosonic creation (annihilation) operator for modes corresponding to long 
wavelength density fluctuations and the operator J has even integer eigenvalues indexing the 
topological winding number of the phase field 4>(z), No = pqL and N is the boson number 



operator. Substituting Eqs. (15) and (16) into (10) yields the mode expanded Hamiltonian 

tcvK 



TXV 



2KL 



J 2 + 



2L 



-(N-N f + J2v\<l\b% 



(17) 



where we have dropped a non-universal constant and defined a new velocity v = ^vjVn 
measuring the famous linear dispersion of the low energy density modes. We have introduced 
the Luttinger parameter K = \Jvn/v^ and although it is well known that no quantum 
phase transition can occur as a function of interactions in a truly Id system described by 



Eq. (17), K can be tuned to initiate a T = crossover between a state with algebraic density 



wave order at K = oo to one with quasi- long range superfluid correlations at K = 0. 

The advantages of having the Hamiltonian in this form are unmistakable due to its 
quadratic nature and the ease with which we can compute averages in the oscillator basis. 
For example, one can exactly determine the grand partition function Z = Tr exp[— (Hll ~ 
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pN)/T) in terms of known special functions^ allowing for the straightforward (although 
possibly tedious) calculation of all two body correlation functions and thermodynamic ob- 
servables in terms of the temperature T, system size L, Luttinger velocity v and Luttinger 
parameter K . We will focus on the derivation of a single observable, the density-density 
correlation function (p(z)p(O)). Such density correlations are of great interest, as they are 
easily measured in numerical simulations and are related via a Fourier transform to the 
experimentally measurable structure factor. In addition, their form provides an intuitive 
qualitative picture of the types of fluctuations (phase or density) which are dominant. This 
knowledge can help to pinpoint which region of the Id crossover phase diagram a given 



system resides in. Starting from the definition of the density operator in Eq. (11) we have 
(keeping only the slowest decaying terms) 

(p(z)p(O)) *p 2 + - 2 (d z 6(z)dA0)) + 2p2(e 2 ^e- 2 ^>. (18) 



Next, using the mode expansion for 9{z) in Eq. (15) the expectation values can be computed 



™ Q -nv/LT 



to glV 



f 2/7 (^-) e - nv / 6LT ] 2/K 
+ Acos(2n Po z)^ d l ^ e _ 7rv/LT) | (19) 

where 6i(x,y) and r](is) are the first Elliptical Theta function and Dedekind Eta function 
respectively. A is a non-universal constant dependent on the short-distance properties of 



the system. Although Eq. (19) may appear daunting at first glance, in the thermodynamic 



limit LT/v — > oo, it simplifies td^ 

(p(z)p(O)) p 2 - + cos(2 W ) (20) 

where it is now clear that K ^> 1 corresponds to a strong tendency towards density wave 
order. 

We are now in a position to test the prediction that for some special radii, the core 
of a helium-4 filled nanopore can be fully described by Luttinger liquid theory. We can 
directly compare the pair correlation function measured in the quantum Monte Carlo with 



the Luttinger liquid prediction of Eq. (19). Before doing so, let us make a brief digression to 



study the connection between the interactions in the truly one dimensional Hamiltonian in 
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FIG. 5. The effective interaction energy between two helium atoms in strictly one dimension 
(lowest curve) as well as for the narrowest (R = 2.9 A) and widest (R = 12.0 A) pores considered 
in this study. The radial density used in the numerical evaluation of Eq. ( |21[ ) was measured at 
T = 0.5 K, but the form of the potential is only weakly dependent on temperature when T < 1.0 K. 

Eq. (j8|) and those experienced by helium atoms inside the quasi-one-dimensional environment 
inside the nanopore. 

B. Effective Interactions and Linear Density in the Inner Cylinder 

The full interaction and external potentials inside the nanopore denned in Eqs. ^ and 
Q and displayed in Figure [7] can be related to the Id potential of Eq. ^ through 

v ld (z) = j J d 2 r J d 2 r'v(r - r')p R (r)p R (r') (21) 

where r = (r, ip, z) is a vector in cylindrical coordinates and Pr(t) is the radial density 
denned in Eq. ([7]). The form of this effective interaction potential can be investigated 
using the measured data from our QMC simulations with the results shown in Figure |5j 
We observe that the effective potential in the R = 2.9 A pore very closely reproduces the 
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FIG. 6. The number of particles per angstrom for helium atoms in strictly one dimension (with 
/j, = 85 K), a pore with R = 2.9 A and inside the inner cylinder with r < ric = 1-75 A for 
R = 12.0 A. All pores have length L = 100 A and the finite radius pores are held at saturated 
vapor pressure with \i = —7.2 K. 

interactions that would be experienced in a system of helium atoms interacting in strictly 
one dimension (R = 0). For the largest radius pore of R = 12.0 A however, screening from 
the concentric shells leads to a potential with a much weaker minimum (about half the depth 
of the unscreened bulk Aziz potential) with its location shifted to smaller separations. The 
effective Id interactions in the axial direction experienced by atoms in the center of the pore 
influences the average separation between atoms which can be determined from the average 
linear density in Figure [6j We observe almost no temperature dependence for small radii and 
the agreement of the linear density for R = and R = 2.9 A is by design, with the chemical 
potential having been tuned to /isvp — ^(0; 2.9) ~ 85 K for the one dimensional (R = 0) 
system to reproduce the energetic confinement displayed in Figure [I] at low temperature. 
The numerical value of po, denned as the zero temperature thermodynamic limit of the 
linear density in the core is measured to be po,o — Po,2.9 = 0.3100(1)A 1 (with the number in 
brackets indicating the uncertainty in the final digit and a second subscript being the pore 
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radius in angstroms). The fact that this value is slightly smaller than r min can be attributed 
to quantum effects. For R = 12.0 A, the observed density of po,i2 = 0.2144(5) A is only 
weakly related to the minimum of the effective Id potential shown in Figure [5] as there is 
a temperature dependent amount of exchange of helium atoms between the inner cylinder 
and surrounding shells. 



C. Density Correlations 

Having studied the effective one dimensional potential felt by helium atoms inside the 
pore, and determined the linear density in the inner cylinder we are now in a position 



to evaluate the efficacy of Eq. (19) in describing density correlations in the nanopore. A 
first glance at this expression indicates that there are four fitting parameters, (po,v, K,A) 
to be determined; a number that would seemingly allow large flexibility (and thus limited 
accuracy) of any least squares fitting procedure. However, one of the most attractive features 



of Eq. (19) is that it provides a prediction for the full analytical form of the finite size and 
temperature scaling behavior of (p(z)p(O)), a rare luxury indeed. We thus have a stringent 
procedure for confirming Luttinger liquid behavior in the nanopore: (1) perform simulations 
for different pore radii, lengths and temperatures measuring thermodynamic estimators and 
correlation functions. (2) The average linear density in the central region of the pore can 
be determined as a function of temperature and extrapolated for each L and R to T = 0, 
this fixes po t R. (3) For each radius which displays a finite density of helium atoms at r = 0, 
perform a least squares fit of the density-density correlation function measured in the QMC 



to Eq. (19) at a single L and T = 0.5 K (the lowest temperature measured). This fixes 
the remaining three parameters, Vr, Kr and Ar where the subscript R will be used to 
distinguish results for different finite radius pores. These parameters are intrinsic to the 



zero temperature thermodynamic limit of Eq. (10) and provided that helium-4 inside the 



core is behaving as a LL, cannot exhibit any temperature dependence. (4) In other words, 



for a given radius, a single fit to Eq. ( 19 ) at fixed temperature and system size is enough to 
compute the LL prediction for the pair correlation function at all other temperatures and 
sizes and we explore the resulting predictions for R = 0, 2.9 and 12.0 A below. Obviously we 
are constrained by the low energy, long wavelength region of applicability of i?LL assumed 
throughout this study and it is easily confirmed that deviations appear at short lengths and 
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FIG. 7. Simulation data (circles) and a fit to Eq. (19) (lines) for the pair correlation function 
of strictly one dimensional helium atoms (R = 0). The chemical potential (fx = 85 K) has been 
set to approximate the energetically confining environment coming from the pore wall (see text). 
Error bars are smaller than the symbol size and data in the main panel has been given a vertical 
T-dependent shift for clarity with T = 0.5 K at the bottom and T = 2.0 K at the top. 

high temperature. 

To orient ourselves we begin by studying the strictly one dimensional system, which we 
fully expect to be well described by LL theory. In Figure [7] we plot the results of QMC 
simulations for the density-density or pair correlation function (symbols) for a chain of 
helium atoms with L = 100 A and /i = 85 K. 

Performing the aforementioned fits at T = 0.5 K, we find vq = 75(1) AK and Kq = 6.4(3) 
where stochastic errors in QMC data produce uncertainty in our regression procedure. The 
exact value of A, is not relevant in the subsequent discussion and it will not be mentioned 
further. Fixing all parameters, Eq. (19) was used to produce the solid lines in Figure (J7|) 
where a vertical shift has been included to allow for the visual differentiation of the effects 
of temperature. We observe spectacular agreement over the entire range of temperatures 
considered, noting that no additional fitting was required above T = 0.5 K. We now perform 
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FIG. 8. Simulation data (symbols) and a fit to Eq. (19) (lines) for the pair correlation function 
along the axis of a nanopore with R = 2.9 A and L = 100 A. Different curves correspond to 
increasing temperature from T = 0.5 K (bottom) to T = 2.0 K (top) (vertical shifts have been 
added for clarity) with the legend displayed in Figure 7. 

an identical procedure for the narrowest finite radius pore, R = 2.9 A, with the results for 
the density-density correlation function shown in Figure [8] Similar to the case of Id helium 
(R = 0), we observe persistent oscillations out to the largest distances possible in the 
L = 100 A pore with periodic boundary conditions indicating a strong tendency towards 
density wave order. As the helium-pore interaction is independent of the axial coordinate z 
and the radial density for R = 2.9 A (Figure [4| exhibits only a single central chain of atoms, 
Galilean invariance further restricts the ratio v 2.9/^2.9 = 7rpo,2.9A^- Extracting values in 
the presence of this constraint at T = 0.5 K yields f 2 .g = 70(3) AK and K 2 .g = 6.0(2), 
in close agreement with the Id chain of helium atoms as expected. For the finite radius 
pore however, we begin to observe deviations from LL predictions both at small distances 
and high temperature. Once the relevant thermal lengthscale, It ~ v/T is on the order 
of the pore diameter, the system can no longer be thought of as quasi-ld and we expect 
significant corrections to arise from the thermal excitation of transverse modes. The values 
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FIG. 9. Simulation data (symbols) and a fit to Eq. (19) (lines) for the pair correlation function 
along the axis of a nanopore with R = 2.9 A for various pore lengths L (increasing from bottom to 
top). As in previous figures, a vertical shift has been added to distinguish the curves and accentuate 
the scaling behavior. Note the use of a dimensionless abscissa in order to plot spatial correlations 
for different pore lengths on the same scale. 

of the LL velocity and interaction strength determined in this way can also be used to test 



the predicted finite size scaling of Eq. (19) at fixed T = 1.0 K as seen in Figure |9j We stress 
that the nearly perfect agreement seen between simulation data and LL theory in Figure [9] 
for L > 75 does not require any additional fitting parameters and is a direct consequence of 
the predictive power of the universal hydrodynamics of Hll- 

Shifting attention to the R = 12.0 A pore, in order to make an effective comparison 
we will focus the analysis on only those helium atoms which spontaneously find themselves 
inside the inner cylinder. The precise definition of which helium atoms are inside the IC is 
somewhat arbitrary, but we define this to coincide with the location of the first minimum 
in the radial density in the presence of finite density at r = 0. For the pores considered 
here, this corresponds to ric = 1.75 A. We expect that the presence of particle exchanges 
between the central chain of helium atoms and the surrounding cylindrical shells should 
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FIG. 10. Simulation data (symbols) and a fit to Eq. (19) (lines) for the inner cylinder pair 
correlation function along the axis of a nanopore with L = 100 A and R = 12.0 A for helium atoms 
with r < nc = 1-75 A. The curves have been shifted for clarity and correspond to increasing 
temperature from T = 0.5 K (bottom) to top T = 2.0 K (top) with with the legend displayed in 
Figure 7. 

begin to play an important role in the 'melting' of the quasi- long range density wave order 
observed in narrow pores. The resulting axial pair correlation function measured in our 



QMC simulations is shown in Figure 10 The data show density-density correlations in the 
inner cylinder that decay much faster than those seen for R = 2.9 A. The lack of Galilean 
invariance due to atoms in the surrounding concentric shells adds some additional freedom to 
the now familiar fitting procedure (and thus more uncertainty) and it has been determined 
that V12 = 42(2) AK and Kyi = 1.3(2). This value of the Luttinger parameters is nearly a 
factor of five smaller than that found for the strictly one dimensional system and indicates 
that this nanopore is in a region of phase space starting to be dominated by superfruid 
fluctuations. 



The success of the Luttinger liquid prediction for (p(z)p(0)) in Eq. (19) in describing the 



density-density correlations measured in our quantum Monte Carlo simulations is convincing 
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TABLE I. Values for the Luttinger velocity and parameter for helium atoms in Id and in the quasi- 
ld core of nanopores by comparing the Luttinger liquid prediction with quantum Monte Carlo 
measurements of the pair correlation function. The number in brackets indicates the uncertainty 
in the final digit. 

Radius R [A] Luttinger Velocity v [AK] Luttinger Parameter K 

0.0 75(1) 6.4(3) 

2.9 70(3) 6.0(2) 

12.0 42(2) 1.3(1) 



evidence that a quantum fluid of helium-4 in pores with R < 3 A can be described by LL 
theory and that the central region of wider pores may display emergent harmonic fluid-like 
properties. In the next section, we discuss the implications of these results and argue that 
the determination of K and v for a given system is not simply an academic exercise but has 
real consequences for the stability of the Luttinger liquid phase. 



V. DISCUSSION 

Although it appears that LL theory works exceedingly well in describing pair or density- 
density correlations for the finite radius pores considered in this study, it is natural to ask 
whether there is something special about this particular observable. As mentioned in Sec- 
tion IV A| the utility of the harmonic fluid description of Id systems is that the resulting mode 



expanded Hamiltonian of Eq. (17) is quadratic (neglecting formally irrelevant operators) and 
thus the full grand partition function can be computed in closed form. We can therefore 
compare other two body correlation functions computed within Luttinger liquid theory with 
those measured in the QMC using the values of Vr and Kr determined in this study and 
displayed in Table [ij We have measured the axial one body matrix n(z) = (^(z)^(0)) in 
the QMCP^l and find acceptable agreement with LL theory (with no new fitting parameters) 
at low T. The large value of K rj 6 for R = and R = 2.9 A adds some complications when 
comparing numerical results for some scalar (non-correlation function) observables with pre- 
dictions coming from bosonization. For example, the strong tendency towards density wave 
order d z 6(z) ~ displayed in Figures [7]{8] points to a compressibility an ~ (iV 2 ) — (N) 2 « 
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and thus an analysis of the probability distribution function for particle number fluctuations 
like the one performed in Reference HH is not feasible. The conjugate relationship between 



density and phase variables in the LL theory (Eq. (12)) would predict that the boson phase 
4>(z) should be totally disordered resulting in a vanishing superfluid density which is observed 
in the numerics. 

For the R = 12 A pore, the analysis presented in the previous section is based on a 
fraction of the total number of 4 He atoms in the pore, those that dynamically have r < ric 
and find themselves in the core of the nanopore. The justification for this originates in 
the idea that we may be able to regard the pore as a coupled multi-component LL, with 
cylindrical shells replacing the legs of previously studied ladders^. Guided by these results 
we assume that only a single gapless degree of freedom may survive as a "center of mass 
mode" in the low energy effective field theory, due to tunneling between the shells. In the 
nanopore, this tunneling has two origins corresponding to the physical mobility of particles 
between shells as well as multi-particle quantum exchange cycles which may dynamically 
connect them at short time scales. 

With this in mind, let us re-analyze the slight discrepancies between the quantum Monte 



Carlo data in Figure 10 and LL theory. At the lowest temperature considered, T = 0.5 K, the 
simulation data appears to show oscillations with a period that is slightly larger than that 
predicted by 1/ (27rpo,i2)- This is most likely attributed to the physical exchange of particles 
between the inner cylinder and the surrounding shells producing a greater uncertainty in 
Po,i2 than is reflected in errorbars and making the extrapolation to zero temperature a 
difficult task. The amplitude of the oscillations in the pair correlation function is also 
slightly overestimated by LL theory at low temperature. Again, the finite radius of the 
pore is to blame, resulting in a finite superfluid fraction of helium (as measured via the 
usual winding number estimator^ along the axis of the pore) of p s /po ~ 0.2 at T = 1.0 K. 
Transverse degrees of freedom that are not frozen out at this temperature lead to an increase 
in multi-particle exchanges that enhance superfluidity in the nanopore. 

The opposite behavior appears to occur at high temperature with the simulation data 
showing weaker decay than predicted by LL theory with thermal fluctuations being unable to 
fully quench the proclivity towards density wave order. This is the opposite effect observed 
for R = 2.9 A in Figure [8] and it can possibly be explained through stabilization of the IC 
from attraction with surrounding atoms. 
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In a previous woriPD finite size corrections to po,i? arising from the inclusion of higher or- 



der formally irrelevant terms in the effective Hamiltonian of Eq. ( 10 ) have been discussed at 
great length. However, due to the relatively large energy scales at play in the nanopore (Fig- 
ure^ they are less important here and are not major players in any observed discrepancies 
between simulation data and the effective field theory. 

We observe a general trend of K decreasing with increasing pore radius as seen in Table [1} 
However, the actual numerical values of the Luttinger parameter K for pores of varying 
radius can provide important information on the sensitivity of the LL to perturbations 
coming from commensuration effects or disorder; both of which are surely present in the real 
experiments of Savard et a/J 14 * 15 l For a strictly one dimensional system, the introduction of 
a weak periodic substrate that is commensurate with the density will only lead to complete 
localization and the destruction of the harmonic fluid if K > 1/2^. Commensuration at 
other wavevectors is less relevant and would require a greater value of K to destabilize the 
LL. The introduction of a weak disorder potential, as might be present near the glassy walls 
of the pore, is known to be relevant only when K > 2/3PI 

For the narrowest pores considered in this study R < 3 A, we have found a value of 
K ~ 6 at saturated vapor pressure, indicating a strong tendency to form a solid, resulting 
from strong confinement and the deep minimum and accompanying hardcore found in the 
effective interaction potential via{z). The experimental confirmation of this result could be 
accomplished by noting that the formation of a quasi-solid should impede superfluid flow 
through a helium-4 filled nanopore at low temperature. 

It may be useful to compare the values of Kr found here with other studies of low 
dimensional helium. For example, in a WA study of helium-4 confined to flow in the 
channels formed by screw dislocations with R ~ 3 A in solid helium, it was found that 
K = 0.205(20)331, nearly thirty times smaller than the comparable value of K 2 .g measured 
in this study. The sources of the discrepancy are rooted in the "softness" of the confin- 
ing potential inside the screw dislocation as helium atoms are able to penetrate into the 
surrounding solid held at /i = 0.02 K corresponding to the bulk melting point. 

Much exciting work remains to be done in the nanopore system including a more sys- 
tematic study of the superfluid density which can be measured in bundles of tubes or pores 
via torsional oscillator technique d 45 * 46 l The construction of more realistic models that con- 
tain both commensuration and disorder potentials would also enhance the applicability of 
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numerical simulations. Further exploration of the available parameter space is also in order, 
including altering the chemical potential at fixed radius to simulate the effects of pressure 
which can be freely tuned in experiments. 

In conclusion, we have studied a quantum fluid of bosonic helium-4 confined inside 
nanopores of varying radii via large scale continuum Worm Algorithm quantum Monte Carlo 
Simulations at saturated vapor pressure below the bulk superfluid transition temperature. 
The results show a progression of phases inside the pore exhibiting a possible quasi-one- 
dimensional core surrounded by concentric shells of helium depending on the radius. When 
the core of the nanopore has a non-zero density of helium, the finite temperature and scal- 
ing properties of the density-density correlation function are fully described within harmonic 
Luttinger liquid theory. The description of the emergent low energy phase of confined helium 
in terms of the harmonic field theory allows for the extraction of the Luttinger parameter K 
which is found to be a decreasing function of radius. As the pore radius increases, the inner 
helium core is screened from the confining effects of the pore wall by the surrounding matter, 
resulting in a pronounced enhancement of quantum exchanges leading to superfluidity. The 
precise relationship between the material through which the pore has been sculpted and the 
resulting Luttinger liquid parameters could be further explored leading to new predictions 
and optimized experiments with the maximum likelihood of detecting a universal and stable 
one-dimensional quantum harmonic fluid of helium at low temperature. 
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